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ABSTRACT 


This manual outlines some of the theoretical aspects embedded in the computer 
code, HITCAN (High Temperature Composite ANalyzer) . HITCAN is a general purpose 
computer program for predicting nonlinear global structural and local stress- 
strain response of arbitrarily oriented, multilayered high temperature metal 
matrix composite structures. This code combines composite mechanics and laminate 
theory with an internal data base of the constituents (matrix, fiber and 
interface) material properties. The thermal and mechanical properties of the 
constituents are considered to be nonlinearly dependent on several parameters 
including temperature, stress and stress rate. The computational procedure for 
the analysis of the composite structure uses the finite element method. HITCAN 
consists of three modules: COBSTRAN, METCAN and MHOST . COBSTRAN generates the 
geometry (pre-processor) and defines the lay up of the different plies. METCAN 
computes the material behavior of the composite and of the constituents. Finally, 
MHOST is a finite element program based on the mixed iterative solution 
technique. MHOST has a library for 2D and 3D isoparametric elements. HITCAN is 
written in FORTRAN 77 computer language and has been configured and executed on 
the NASA Lewis Research Center CRAY XMP and YMP computers. 
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CHAPTER 1 
INTRODUCTION 


HITCAN (High Temperature Composite ANalyzer) is a general purpose computer program 
for predicting the global structural response of high temperature metal matrix composite 
structures. This document describes the basic theory behind this code. 

The HITCAN computer program is a combination of three computer codes developed in-house 
at the NASA Lewis Research Center. They are: 

- a mesh generator; COBSTRAN (Reference 1); 

- a metal matrix composite analyzer; METCAN 
(Reference 2); 

- a finite element structural analysis code; MHOST 
(Reference 3). 

In Chapter 2, the theoretical aspects of the METCAN-MHOST synthesis will be 
discussed. The capabilities and limitations of HITCAN are discussed in Chapter 3. In 
Appendices 1 and 2, the theory behind the individual library modules METCAN and MHOST 
will be reviewed. Since the functionality of COBSTRAN as used in HITCAN, is restricted to 
only the mesh generation, its theoretical aspects will not be discussed here. The reader is 
referred to Reference 4 for further information. 
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CHAPTER 2 

METCAN-MHOST SYNTHESIS 

As mentioned previously, HITCAN is a computer program for predicting the global 
structural response of high temperature metal matrix composite structures. At high temperature, 
metal matrix composite structures experience nonlinear material behavior. As a result of this 
nonlinear behavior, the finite element equations are nonlinear. Consequently, for HITCAN to 
be of any use it must be able to solve these equations. Thus, a nonlinear solution scheme was 
incorporated into HITCAN. 

The nonlinear solution scheme consists of two parts. The first part is an analysis of the 
composite material at the laminate, ply, and constituent levels. To perform this local analysis, 
METCAN was chosen. The second part of the scheme is a global analysis of the structure using 
MHOST. Thus, the core of this nonlinear solution scheme is the METCAN-MHOST synthesis. 
Both METCAN and MHOST are stand-alone computer programs. Hence, it was desirable that 
neither one of these programs would be modified extensively in order to work in HITCAN. In 
order to minimize the number of modifications, it was decided that the finite element formulation 
for linear structural analysis would serve as the basis for nonlinear structural analysis. For the 
nonlinear analysis, the linear formulation is applied in an incremental/iterative manner. Before 
proceeding further, a brief review of the finite element formulation for linear analysis is in 
order. 


The global formulation of a linear static structural analysis problem using the finite 
element method is given in the following expression: 

[K]{u} = {F} 
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In this expression, [K] is the global structural stiffness matrix, {u} is the vector of nodal 
displacements, and {F} is the vector containing the applied nodal loads. 

The next expression is the eigenvalue problem corresponding to a free vibration analysis. 

([K] - *[M]){u} = {0} - {6> n } 

In this expression, [K] is the global structural stiffness matrix, [M] is the global structural 
mass matrix, k is the eigenvalue associated with the solution of the problem, [I] is the identity 
matrix, and {o n } is vector containing the natural frequencies where n is the number of modes. 

The last expression is the eigenvalue problem corresponding to a static buckling analysis. 

([K] + «KJ){u} = {0} - {S CT }, 

where [K„] is the stress stiffness matrix, [K] is the linear stiffness matrix, k some scaler 
multiplier of [KJ, {0} is the null column vector, and {S CT } is the critical buckling load. 

Although the formulation summarized above is for a general linear structural analysis 
problem, it can be applied in the nonlinear analysis of a structure. How this can be 
accomplished is described in the following section. 
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2.1 NONLINEAR ANALYSIS 

The finite element formulation for a linear structural analysis can serve as the basis for 
a nonlinear structural analysis. For a nonlinear analysis, the linear formulation is applied in a 
incremental/iterative manner. This approach allows for the fact that the material and thermal 
properties are not constant over the range of the problem, but are actually dependent functions 
of the primitive variables, such as stress, rate of stress, temperature, rate of temperature, 
metallurgical reaction between fiber and matrix, and time. 

In applying this incremental/iterative approach, the solution of a nonlinear structural 
analysis problem is treated as an iterative process that involves solving the linear problem 
formulation for distinct load increments of the total range of the solution for the problem. 
Furthermore, as the iterative solution process progresses, the results from the analysis of each 
step or increment are used to update the material and thermal properties of the composite at each 
successive increment. The flow diagram for the incremental/iterative approach is shown in 
Figure 2.1. 

In HITCAN, the iterative procedure chosen was the direct iteration method. Assuming 
a single degree-of-freedom system and that the stiffness is a function of the temperature, 
temperature rate, stress, and stress rate, the equation that needs to be solved is [K(T,T, or, cr, t)JD 
= F. The direct iteration method works as follows. 

- Initially, the material properties are calculated based on the nodal temperature, 
temperature rate, stresses, and stress rates obtained from the previous load increment 
and time. 
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- The incremental displacements are calculated from the equation 
AD 1 = [K(T ,T, a, <r ,t)°]‘ ! AF, 

where [K(T,T,cr,o',t)°] is the initial slope and AF is the incremental load. 

Repetition of this process can be written as 
AD" = [KCTJ^^^^TAF. 

Here the slope [K(T,T,ff,cr,t) n ' 1 ] is calculated based on the material properties of the 
previous iteration. The material properties are functions of the nodal temperature, temperature 
rate, stresses, and stress rate. 

- This process is terminated when the error becomes 
sufficiently small, i. e. , 

1 1 AD n - AD”" 1 1 1 / 1 1 AD 11 " 1 - AD 11 " 2 1 1 < tolerance. 

This process is shown graphically in Figure 2.2. 
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F 



DISPLACEMENT, D 


DIRECT ITERATION METHOD 
FIGURE 2.2 
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2.2 COMPOSITE ANALYSIS 

The distinguishing features involved in the analysis of a composite structure include the 
additional theories and assumptions employed to represent the inherently heterogeneous 
composite structure as an equivalent "pseudo-homogeneous” structure. In the approach taken 
in HITCAN, these additional theories provide for the generation of the various global structural 
properties and for the degeneration of the material properties during the incremental solution 
process. These aspects, unique to the analysis of a composite structure, involve such subjects 
as composite micromechanics theory and laminated plate theory. 

The iterative approach used in HITCAN, as illustrated in Figures 2.1 and 2.2, is shown 
schematically in Figure 2.3. As can be seen in the figure, the direct iteration method has been 
modified to account for the heterogeneity of composite materials. The key elements comprising 
this analysis include: 

multifactor interaction model for nonlinear material; 

- composite micromechanics theory; 

laminate theory; 

finite element structural analysis. 

The circular arrangement of Figure 2.3 is intended to illustrate the iterative solution 
strategy imbedded in the incremental/iterative approach as it is used in HITCAN. Figure 2.4 
shows the quantities passed between METCAN and MHOST, namely, stress resultants, laminate 
properties, thermal strains, and composite properties. 
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As can be seen in the Figure 2.3, the two library modules METCAN and MHOST 
comprise the iterative process. In this scheme, iteration is required at two levels. The first or 
local level is in METCAN. This is required since, in the nonlinear material model, the 
constituent properties are functions of the constituent temperature, temperature rates, stresses, 
stress rates, and time. The second or global iteration is required to establish equilibrium of the 
global structural response. 

From Figure 2.3, it can be seen that METCAN performs steps E, F, A, B, and C; i. e., 

- METCAN takes the stress resultants obtained from the structural response and 
determines the ply stresses, using laminate theory; 

- Using composite micromechanics theory METCAN determines the constituent stresses 
and stress rates; 

- METCAN then calculates new constituent material properties using the multi-factor 
model; 

- From micromechanics theory, METCAN is then able to calculate new ply properties; 


- With the new ply properties, METCAN then determines the new laminate properties; 


- Using laminate theory, METCAN calculates the strains at each node. Using these ply 
strains, convergence is checked within METCAN. 
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The iterative process occurring in METCAN comprises the steps E, F, A, B, and C. Since 
these steps are contained entirely within METCAN, it is called only once for each node. As 
shown in Figure 2.3, MHOST simply determines the global response. If a convergence criteria 
is included in the loop, Figure 2.3 can be redrawn as shown in Figure 2.5. Figure 2.6 shows 
the flow diagram of the incremental/iterative process used in HITCAN (subroutine NLINER). 
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CHAPTER 3 

CAPABILITIES AND LIMITATIONS 


HITCAN is capable of predicting global structural and local stress-strain response of 
multilayered HTMMC structures exhibiting nonlinear material behavior. The constitutive model 
employed in METCAN is specifically designed for HTMMC, therefore it is recommended that 
HITCAN be used only for metal matrix composites. 

At the present time, the following analyses are available in HITCAN: 

- Incremental static analysis with nonlinear anisotropic material behavior; 

- Dynamic analysis using direct time integration; 

- Modal analysis (free vibration); 

- Buckling analysis (first critical buckling load). 

The element library includes 3 four-noded elements, i.e. plate, plane stress, and plane 
strain, and 1 eight noded element, i.e. a 3D solid. The current mesh generation capability of 
the code allows modeling of solid structures using any of the 4 types of elements. The user may 
also input a finite element model directly. The code is capable of handling a variety of boundary 
conditions, loadings (centrifugal, concentrated, distributed, pressure, temperature, static, 
transient, cyclic, and impact), and various types of structures (such as beam, plate, ring, curved 
panel, and built-up structures). A list of HITCAN’ s capabilities can be found in Table 3.1. 
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The limitations of the code are : 

- Formulation assume small displacement and small strain theory; 

- Elements of different types cannot be combined; 

- Hollow structures can be modeled using the plate element only; 

- The finite element model generated by HITCAN can have a nonuniform mesh only 
along the x-axis; 

- If the curvature is large, the mesh will not be uniform. 

HITCAN presents analysis results at the global, laminate, and ply levels. Results include 
displacements, reactions, stresses, and modes shapes. The code also has the capability to 
generate post-processing files for PATRAN. 
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““ . Type of Element -«► 

Type of Analysis^ ■ — 

Plate 

Plane 

Strain 

Plane 

Stress 

8-Node 

Solid 

Static 

tested 

- 

- 

tested 

Buckling (a) 

tested 

- 

- 

- 

Load Stepping 

tested 

- 

- 

tested 

Modal (Natural Vibration Modes) (b) 

tested 

- 

- 

- 

Time-domain 

- 

- 

- 

- 

Loading 

Mechanical 

tested 



tested 

Thermal 

tested 



tested 

Cyclic 

- 

- 

- 

- 

Impact 

- 

- 

- 

- 

Constitutive Models (c) 
P = Constant 

tested 



tested 

P = f( T ) (temperature dependence) 

tested 

. 

- 

tested 

P = f( a ) (stress dependence) 

tested 

- 

- 

tested 

P = f (or) (stress rate dependence) 

tested 



tested 

p - f ( t ) (creep) 

- 

. 


- 

P = UTd , a a) (combination) 

tested 

_ 

_ 

tested 

p = f( Tc, <j ,a , t ) (creep combination) 

- 


_ 

- 

Fiber Degradation 

tested 


_ 

tested 

Fabrication-induced Stresses 

tested 



tested 

Ply Orientations 
Arbitrary 

tested 



tested 


(a) Tested 1 buckling mode (c) Constitutive models: Notation 
<b) Tested 4 vibration modes £• Material properties a : Stress 

T: Temperature or : Stress rate 

t : Time 


TABLE 3.1: HITCAN Capabilities for Composite Materials 
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APPENDIX 1 
METCAN THEORY 

In this chapter, the theoretical aspects of METCAN will be reviewed. Only a brief review 
can be given here. For more detailed information the reader is directed toward References 5 - 
10 . 


As stated in the previous chapter, METCAN performs the following iterative functions in 
HITCAN: 

METCAN takes the stress resultants obtained from the structural response and 
determines the ply stresses, using laminate theory; 

using composite micromechanics theory METCAN determines the constituent 
stresses; 

METCAN then calculates new constituent material properties using the multi- 
factor model; 

from micromechanics theory, METCAN is then able to calculate new ply 
properties; 

with the new ply properties, METCAN then determines the new laminate 
properties; 

using laminate theory METCAN, calculates the ply strains at each node. These 
strains are used in checking convergence. 
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It was also stated previously that because of the form of the nonlinear constitutive model, 
iterating is required within METCAN. With all of the above in mind, the flow diagram for 
METCAN can be drawn as shown in Figure Al.l. For the check on convergence, METCAN 
uses the ply strains; i. e. , 

(e?" - e i ° Id )/€ i okl < tolerance, 
where i = 1,...,6 at each ply. 

The rest of this chapter is devoted to the key elements of METCAN, namely, composite 
micromechanics theory, laminate theory, nonlinear constitutive model, and fiber degradation. 
Composite micromechanics theory can be found in Appendix 1.1. In Appendix 1.2, a brief 
review of laminate theory is given and in Appendix 1.3 the constitutive model and fiber 
degradation are discussed. 
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METCAN FLOW DIAGRAM 
FIGURE A1.1 
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APPENDIX 1.1 - COMPOSITE MICROMECHANICS THEORY 

Composite micromechanics theory relies on the principles of solid mechanics, 
thermodynamics, etc. at different levels of mathematical sophistication, together with certain 
assumptions and approximations. In the approach taken in METCAN, the principles of 
displacement compatibility, force equilibrium, and Fourier’s law for heat conduction are used. 
In addition, the following assumptions are made: 

- fibers are continuous and parallel; 

- properties of all fibers are identical; 

- complete bonding exists between constituents. 

The periodic structure of an unidirectional metal matrix ply is approximated by a square 
unit cell model. The geometry of the model is shown in Figure A 1.2. It should be noted that 
the interphase growth is assumed to result from the degradation of fiber material and thus 
propagates inward causing a continuous decrease of the current fiber diameter from the original 
fiber diameter. With the existence of the interphase, three subregions (A, B, C) are 
distinguished to characterize the intralaminar (through-the-thickness) nonuniformity of the 
constituent microstresses and material properties. Further information on the composite 
micromechanics can be found in Reference 5. 
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APPENDIX 1.2 - LAMINATE THEORY 

In this section, a brief review of laminate theory will be given. For a more complete 
discussion of laminate theory, the reader is referred to Reference 12. 

Laminate theory provides a formal procedure to relate the behavior of a laminated 
composite structure to the behavior of the individual laminate of the composite. An 
unidirectional fiber-reinforced lamina can be approximated as a homogeneous orthotropic ply. 
From assumptions analogous to the Kirchoff hypothesis in classical plate theory, the laminae of 
a laminated composite structure are treated as thin plates and therefore are assumed to exist in 
a plane stress state. As such, the constitutive relationship for the k* specially orthotropic ply 
(referred to the lamina material coordinate system as shown in Figure A 1.3) is given by 


'I 

6 11 


E 11 E 12 0 

-1 

( \ 
a u 

< c 22 

V = 

E 12 E 22 0 


°22 ” 

, e 12. 

k 

KO 

W 

o 

o 

1 

k 

^12, 


In terms of stress, this can be rewritten as 
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Appendix 1.2 March, 1992 




LAM/THEORY3 



HLTCAN Theoretical Manual - Version 1.0 

where the subscript m indicates ply material coordinate system. For the generally orthotropic 
ply (referred to in the global structural coordinate system) the constitutive relationship is given 
by; 



r \ 

r \ 


a X2 




a 

YY 

> * ^l£ a k ^2m^k " 

E yy y 



k 

Z XV 

\ * J 


where the subscript s indicates global structural coordinate system. The matrix [RJ t is the 
rotation or transformation matrix and is defined as; 


(R 



k 


2 

cos 0 

sin 2 - 

%sin29 

si n 2 0 

2 a 
cos 0 

- i isin29 

-sxr.20 

sin29 

cas25 


and the angle 0 represents the orientation of the lamina material coordinate axes (1, 2) for the 
ply with respect to the global structural coordinate axes (x, y). 
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In terms of the midplane strain and laminate curvature, the constitutive relationship for the k* 
generally orthotropic lamina becomes; 


{0 i3 } k = 


[E 


2s J k 


l£ is ; k 1 Zs d k 




0 \ * 

V*C }) 


where the subscript c denotes composite and z°, is the perpendicular distance from the laminate 
midplane to the center of the k* ply. 


The stress resultants (N x , N y , N v , M*, My, M^) are shown in Figure A1.4. The resultants are 
given by the following equations; 


r \ 
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The single continuous integrals above can be separated into sums of integrals over each ply. 
The stress resultants are then given by; 




n 

E 

k=I 



{a 2s'k 


dz 


{M } 
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n 
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k=I 
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z dz 
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These equations can be rewritten as 
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These equations can be simplified to 


where, 


<N C } = £A_; + [C c ] {,c c } 

(•V) = (e°J + [D ] {<„} 


n 


- C CS.-J 


c- ** " Zjc * i) 


[C J = 4 : [e b .] v (z£ - z,L.) 


- k~l 
n 


Zs J k K-l 

3 _3 


lD c ] “ f (z k “ 2 k-l } 


Although the through-the-thickness stresses (a a , cr^, a«) are neglected in deriving the 
constitutive relationships for the orthotropic lamina, these stresses generally will exist in the 
laminate. These stresses are determined from equilibrium considerations and can be written as 
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Assuming a parabolic distribution of shear stress (a^, <?„) over the plate thickness, the shear 
stress for the ply can be written as 


’txz 


1 R* 

4 h 


i- 


tyz 


2 Q ' 
4 


y 

h 


<ir> 


The through-the-thickness normal stress for the ply can be written as 


Izz 


= ( 




j>&) ( 2 r f? 

d y ’ K 4 1 h 
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APPENDIX 1.3 - MATERIAL MODEL AND FIBER DEGRADATION 


A modular multifactor constitutive model accounting for the effect on constituent properties 
of several parameters such as temperature, stress, and stress rate is employed. This model is 
shown below. 


(1) Mechanical property (moduli, strength) P M 



p M 

r t„ - t~i 

n 

S F - a 

m 

s F - 

I 

P Mo 

Jm - T 0 J 


s F - <^o_ 


Sp - (T 


(2) Thermal property (expansion coefficients, thermal conductivity, heat 

capacity) P T 






Pt 
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S F - a 


S F - fro 









where 


M 
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Me > P To 
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T< 

S l 
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a 

n, m, 1 


denotes the current mechanical property of interest 
denotes the current thermal properly of interest 
are corresponding properties at reference conditions 
is the melting temperature 
is the current temperature 

is the reference temperature at which P Me & Ep o are determined 

is the fracture stress determined at T conditions 

o 

is die current stress 

is the reference stress at which P Me & Bp 0 are determined 

is an appropriately selected stress rate, for example, the stress rate at which 
penetration occurs during impact 

is the stress rate at which P Me &E^ o are determined 

is the stress rate, and 

are empirical constants 


Appendix 1.3 


March, 1992 


MAT/FIB/DEGRADl 




HITCAN Theoretical Manual - Version 1.0 


Fiber degradation is a metallurgical phenomenon common in the refractory-fiber metal- 
matrix class of high temperature composites. 

The phenomenon involves a chemical reaction at the fiber-matrix interface of a composite 
lamina. The exact nature and extent of the degradation process is dependant on the particular 
combination of fiber and matrix in the composite system. The result of this degradation process 
is the creation of an interphase of material at the original fiber/matrix boundary. 

Depending on the propensity for a particular fiber/matrix system to interact and the 
character of the resulting interphase of material, the fiber degradation can have a significant 
deteriorating effect on the properties of the composite. Figure A 1.6 illustrates the fiber 
degradation process used in METCAN. In the expression: D 0 represents the initial fiber 
diameter, D the fiber diameter after degradation, and F is the percent of fiber degradation. 
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FIBER DEGRADATION IN METCAN 
FIGURE A1. 6 
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APPENDIX 2 
MHOST THEORY 


As shown in Figure 2.3, the MHOST computer program is used to perform the global 
structural analysis. MHOST is a solid and structural analysis program based on mixed finite 
element technology, and is specifically designed for three-dimensional inelastic analysis. 
MHOST includes two- and three-dimensional elements along with beam and shell structural 
elements. Many options are available in the constitutive equation library, the solution 
algorithms, and the analysis capabilities. The options chosen in MHOST for use in HITCAN 
are: 


- material option composite laminate for use with the plate element; 

- material option anisotropic elasticity for the plane stress, plane strain, and 
eight-node solid elements; 

- linear static, buckling, and modal analysis for the analysis option; 

- displacement option for the finite element formulation. 

In the following sections, brief theoretical reviews are given of those aspects of MHOST 
used in HITCAN. For further information the reader is directed to References 13, 14, and 15. 
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APPENDIX 2.1 - DISPLACEMENT FORMULATION 

Using the displacement based formulation, the finite element problem can be reduced to the 
matrix equation 


[K]{u} = {F}, 


where 


[K] = J 0 [B] T [D][B]dQ 

and where [B] and [D] are, respectively, the strain-displacement operator and the stress-strain 
relation for the element. 

A solution is obtained by solving for the vector {u} of unknown nodal displacements. 
Strains and stresses are not solved for directly. In MHOST, the strains are determined using 

{e c } = [Bj{u c }, 

where {ejis the strain defined at the integration points within the element and {u e } is the 
element nodal displacement vector. These strains are then extrapolated to the nodes. Here the 
stresses are calculated using 

{*} = IP]{e}, 

where {a} is the nodal stress vector, [D] is the stress-strain relation, and {e} is the nodal strain 
vector. 
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APPENDIX 2.2 - ELEMENTS 

The element library includes 3 four-node elements; i. e., plate, plane stress, and plane 
strain, and 1 eight-node element, a 3D solid. This section gives the theoretical background for 
these elements. 

The plane stress element is a bilinear isoparametric element, based on selective reduced 
integration techniques. The element is a four noded quadrilateral with nodes numbered in a 
counter-clockwise fashion as shown in Figure A2.1. The position of each node is defined by 
a {x, y} coordinate pair. 

The displacements, strains, and stresses are stored as a set of vectors with components 
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f °xx 
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*yy 
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°yy 

“rl 


Yzy i 
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where u, and u, are the two components of displacement, e x and e y are the two direct strain 
components, y^ is the shear strain (using the engineering definition), cr x and cr y are the direct stress 
components, and is the in-plane stress. 

The strain-displacement relation at a point x 5 within the element is defined by 

D 

£/=B/u*= X 

a=A 
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PLANE STRESS ELEMENT 
FIGURE A2.1 
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where the {uj for a = A, B, C, and D are the nodal displacement vectors, 
be partitioned into two submatrices 


where [By are the terms associated with the direct strain components 


f^a 



0 



and [B S J is used to compute the shear strain component 


r dN , dN a 


x - x ± 


x J 


The stress-strain relation used is 


[D] = 


^11 •^'llM’12 Q 

^11 1*12 ^22 q 

0 0 G 12 
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The plane strain element is a bilinear isoparametric element. The element is a four-node 
quadrilateral with nodes numbered in a counter-clockwise fashion as shown in Figure A2.2. 

The displacements, strains, an stresses for this element are stored as a set of vectors with 
components 



where u* and Uy are the two components of displacement, e x , e y , and e z are the direct strain 
components (with e z =0), y^ is the shear strain (using the engineering definition), a x , <r y , and a z 
are the direct stress components, and is in-plane shear stress. 

The strain-displacement relation at a point x ; within the element is defined by 

D 

e,= 8,u«= £ B/ a u a 

a =i4 
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PLANE STRESS ELEMENT 
FIGURE A2.2 
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where the {u,} for a ~ A, B, C, and D are the nodal displacement vectors, 
the element can be partitioned into two submatrices 
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ta 
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where [B d J are the terms associated with the direct strain components 


Ml 

Bx 
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1 
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By 


0 0 


X = Xi 


whereas [B‘J is used to compute the shear strain component. 
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The stress-strain relation used is 


^11 ( M > 21 + M , 3lM-23 J ^ll^ 

^1 ( 1 — 1 ^ 23 ^ 32 ^ 


•^22 ( jj- ^31^31 ) 
^”l J '3lM’32^ 
0 
0 


^22 

C 2 

0 

0 


0 

1 

0 



where 


C l = + (l i 12 + l A 13^32)/( 1 -^23i A 32)^ 

^2 = ^ ~M’23l A 32^ t ( 1^21 + ^31^23 ) ^12 + ^13 ^32 ^ ^ ”^13^31^ 
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The plate element is a bilinear isoparametric variable-thickness plate element based on 
the Reissner-Mindlin plate and shell theories. These theories account for shear deformation, and 
can be applied to moderately thick shells. 

The element is a four-noded quadrilateral in three-dimensional space as shown in Figure 
A2.3. Each node is defined by four coordinates {x, y, z, t} where t denotes the thickness of the 
shell at each node. 

The plate element is formulated in terms of global translations and rotations, generalized 
strains, and stress resultants. These displacements, generalized strains and stress resultants are 
stored as a set of vectors with the components 
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PLANE ELEMENT 
FIGURE A2.3 
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where u,, Uy, and u 2 are the three translational displacements, 0 X , 0 y , and 0 2 represent three 
rotations, all expressed in the global Cartesian coordinate system; e x and e y are two in-plane 
strains, Yx y is in-plane shear strain, and are the transverse shear strains, k*, Ky, and k, 
denote the three curvatures at a point, and Y z is the term associated with the "drilling'' rotation, 
all expressed in the local element coordinate system; N x and N y are the in-plane membrane 
forces, the in-plane shear, and are two transverse shears, M*, M y , and are the 
bending moments at a point, and the force resultant associated with the "drilling" mode, all 
expressed in the local element coordinate system. 

As mentioned previously, the plate element is a bilinear isoparametric variable-thickness 
element based on the Reissner-Mindlin plate and shell theories. The weak variational form for 
these equations involves only first-order derivatives of displacement, permitting the use of linear 
interpolation for both rotations and transverse displacements. The strain-displacement relation 
for an arbitrary point x ; at the mid-surface of the plate is defined by 

D 

< 2 , = B lU * S 

a - A 


where the {u a } for a = A, B, C, and D are the nodal displacement vectors. The [B] matrix for 
the element can be partitioned into four submatrices 


Appendix 2.2 


March, 1992 


ELEM10 



HITCAN Theoretical Manual - Version 1.0 


Bia- 


dN a 


0 0 0 0 0 
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where [B m jJ denotes the terms associated with the two direct membrane strain components 
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where [B’J denotes the terms associated with the two direct transverse shear strain components 


Bia = 


BN . 

0 0 0 0 — 0 
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where [B b J denotes the terms associated with the two direct bending strain components 

dN a BN a 






0 0 0 -2 N a 


where [B d J denotes the terms associated with the two direct drilling strain components 
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The stress-strain relation for this element is entered into MHOST using the axial, bending, 
and coupling matrices of laminate theory. 

The eight-node solid element is a trilinear isoparametric element. The element is shown 
in Figure A2.4. The position of each node is defined by a {x, y, z} coordinate vector. 

The displacements, strains, and stresses of this element are stored as a set of vectors with 
components 
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where u*, Uy, and u 2 are the three components of displacement, e x , e y , and e z are the direct strain 
components, y w , y^, and y^ 

are the three shear strains (using the engineering definition), a x , cr y , and a z are the direct stress 
components, and r^, and 
are the three shear stresses. 
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EIGHT-NODE SOLID ELEMENT 
FIGURE A2.4 
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The strain-displacement relation at a point x ; within the element is defined by 

H 

£»=B,u«s £ B faUa 
a=A 


where the {u,} for a = A, B, H are the nodal displacement vectors. The [B] matrix for the 
element can be partitioned into two submatrices 
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where [B'J are the terms associated with the direct strain components 
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whereas (B‘J is used to compute the shear strain components 
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The stress-strain relation used is 
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where C - ^ 12^21 ^ 23^32 ^ 31^13 ^^ 21 ^ 32^13 
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APPENDIX 2.3 - BUCKLING AND MODAL ANALYSIS 

In this section, a brief review of the theoretical aspects of buckling analysis and modal 
analysis is given. In MHOST both of these analyses are formulated as eigenvalue problems. 
The subspace iteration method is used in MHOST to determine the eigenvalues and eigenvectors. 

The eigenvalue problem associated with the critical buckling load can be written as 

([K] + X[KJ){u} - {0}, 

where [KJ is the stress stiffness matrix, [K] is the linear stiffness matrix, and X some scaler 
multiplier of [KJ. Since MHOST, as it is used in HITCAN, only performs a linear analysis, 
then neither [K] or [KJ are functions of displacement. The critical buckling load is associated 
with A. ct , the lowest magnitude eigenvalue. The displacement vector {u} identifies the buckled 
shape, but not its magnitude. For further information, see Reference 12. 

The eigenvalue problem associated with modal analysis is 

m - mm]){u} = {o}, 

X = or, 

where [K] is the stiffness matrix, [M] is the mass matrix, {u} is the amplitude of the nodal 
degrees-of-freedom, o> is the frequency, and X is the eigenvalue. Since neither geometric nor 
material nonlinearity are used in MHOST, then neither [K] or [M] is a function of o>. Only 
certain values X t will satisfy the above 
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equations. These A. ; are called eigenvalues. To each eigenvalue there corresponds an 
eigenvector {u ; }. The eigenvalue problem is to extract the solution pairs A.; and {uj. 

As mentioned previously, the subspace iteration method is used to solve the eigenvalue 
problem. For a description of this method, see Reference 14. 
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